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Counting toroidal binary arrays, II 

S. N. Ethier* and Jiyeon Lee^ 


Abstract 

We derive formulas for (i) the number of toroidal nxn binary arrays, 
allowing rotation of rows and/or columns as well as matrix transposition, 
and (n) the number of toroidal n x n binary arrays, allowing rotation 
and/or reflection of rows and/or columns as well as matrix transposition. 
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1 Introduction 

A previous paper [1] found the number of (distinct) toroidal mxn binary arrays, 
allowing rotation of rows and/or columns, to be 

a(m,n) := - V V <p(c)<p(d) 2 mn / lcm ^ d \ (1) 

mn ' z — 1 ' 

c | m d\n 

where ip is Euler’s phi function and 1cm stands for least common multiple. This 
is A184271 in the On-Line Encyclopedia of Integer Sequences [2]. The main 
diagonal is A179043. It was also shown that, allowing rotation and/or reflection 
of rows and/or columns, the number becomes 

b(m, n ) := 6i(m, n) + 62(771, n) + 63(771, n) + 64(771, n), (2) 

where 

64(771, n) ^ X X^W2 m " /lcm(c ’ d) , 

c | m d\n 
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b 2 (m , n) 



d | n 



2-2 mn / d ), if mis even, 


if m is odd; 


b 3 (m,n) := b 2 (n,m), 


and 


& 4 (m, n) : 


2(mn-3)/2^ m an d n are odd; 

3 • 2 mra / 2-3 , if m and n have opposite parity; 
7 • 2 m ™/ 2-4 , if m and n are even. 


(The formula for b 2 {m , n) given in [1] is simplified here.) This is A222188 in the 
OEIS [2]. The main diagonal is A209251. 

Our aim here is to derive the corresponding formulas when m = n and we 
allow matrix transposition as well. More precisely, we show that the number of 
(distinct) toroidal nxn binary arrays, allowing rotation of rows and/or columns 
as well as matrix transposition, is 



( 3 ) 


where a(n,n ) is from (1). When we allow rotation and/or reflection of rows 
and/or columns as well as matrix transposition, the number becomes 


p(n) = i b(n, n) + -^ J2 v( d ) 2 n( ^ n+d - 2 ^ 2 ^^ 2d '> 


d | n 


2 (n 2 -5)/4 ) if n i s odd; 
5 . 2 n / 4_3 ) if n is even, 


( 4 ) 


where b(n,n) is from (2). At the time of writing, sequences (3) and (4) were 
not in the OEIS. 

For an alternative description, we could define a group action on the set of 
nxn binary arrays, which has 2" elements. If the group is generated by a 
(row rotation) and r (column rotation), then the number of orbits is given by 
a(n,n); see [1]. If the group is generated by a, r, and £ (matrix transposition), 
then the number of orbits is given by a(n); see Theorem 1 below. If the group 
is generated by tr, r, p (row reflection), and 9 (column reflection), then the 
number of orbits is given by 6(n,n); see [1], If the group is generated by er, r, 
p, 9, and </, then the number of orbits is given by /3(n); see Theorem 2 below. 
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Both theorems are proved using Polya’s enumeration theorem (actually, the 
simplified unweighted version; see, e.g., van Lint and Wilson [3, Theorem 37.1, 
p. 524]). 

To help clarify the distinction between the various group actions, we consider 
the case of 3 x 3 binary arrays as in [1], When the group is generated by a and 
r (allowing rotation of rows and/or columns), there are 64 orbits, which were 
listed in [1], When the group is generated by er, r, and £ (allowing rotation 
of rows and/or columns as well as matrix transposition), there are 44 orbits, 
which are listed in Table 1 below. When the group is generated by er, r, p, 
and 9 (allowing rotation and/or reflection of rows and/or columns), there are 
36 orbits, which were listed in [1]. When the group is generated by a, r, p, 9 , 
and </ (allowing rotation and/or reflection of rows and/or columns as well as 
matrix transposition), there are 26 orbits, which are listed in Table 2 below. 


Table 1: A list of the 44 orbits of the group action in which the group generated 
by cr, t, and £ acts on the set of 3 x 3 binary arrays. (Rows and/or columns 
can be rotated and matrices can be transposed.) Each orbit is represented by 
its minimal element in 9-bit binary form. Subscripts indicate orbit size. Bars 
separate different numbers of Is. 


/0 0 0 \ 
0 0 0 

V 1 1 Vs 
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Table 2: A list of the 26 orbits of the group action in which the group generated 
by cr, r, p, 9 , and ( acts on the set of 3 x 3 binary arrays. (Rows and/or columns 
can be rotated and/or reflected and matrices can be transposed.) Each orbit 
is represented by its minimal clement in 9-bit binary form. Subscripts indicate 
orbit size. Bars separate different numbers of Is. 



Table 3 provides numerical values for a(n ) and /3 (n) for small n. 
Table 3: The values of a(n) and j3{n ) for n = 1, 2,..., 12. 


n <*(«) 3(n) 

1 

2 

2 

2 

6 

6 

3 

44 

26 

4 

2209 

805 

5 

674384 

172112 

6 

954623404 

239123150 

7 

5744406453840 

1436120190288 

8 

144115192471496836 

36028817512382026 

9 

14925010120653819583840 

3731252531904348833632 

10 

6338253001142965335834871200 

1584563250300891724601560272 

11 

10985355337065423791175013899922368 

2746338834266358751489231123956672 

12 

77433143050453552587418968170813573149024 

19358285762613388352671214587818634041520 


We take this opportunity to correct a small gap in the proof of Theorem 2 
in [1]. The proof assumed implicitly that to, n > 3. The theorem is correct as 
stated for to, n > 1, so the proof is incomplete if to or n is 1 or 2. Following the 
proof of Theorem 2 below, we supply the missing steps. 
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2 Rotation of rows and columns, and matrix 
transposition 

Let X n := {0, l j{o,i,...,n-i} 2 b e the set ofnxn matrices of Os and Is, which has 
2 n elements. Let a(n) denote the number of orbits of the group action on X n 
by the group of order 2 n 2 generated by a (row rotation), r (column rotation), 
and £ (matrix transposition). (Exception: If n = 1, the group is of order 1.) 

Informally, a(n) is the number of (distinct) toroidal n x n binary arrays, 
allowing rotation of rows and/or columns as well as matrix transposition. 

Theorem 1. With a(n,n ) defined using (1), a(n) is given by (3). 

Proof. Let us assume that n > 2. By Polya’s enumeration theorem, 

1 n— 1 n —1 

= ( 5 ) 

2=0 j—0 

where Aij (resp., Eij) is the number of cycles in the permutation a l r J (resp., 
crV^C); here a rotates the rows (row 0 becomes row 1, row 1 becomes row 2, ..., 
row n — 1 becomes row 0), r rotates the columns, and £ transposes the matrix. 
We know from [1] that 


a(n,n) = -^^^2 A «, (6) 

2=0 j= 0 

so it remains to find E^. The permutation </ has n fixed points and (”) trans¬ 
positions, so E 00 = n{n + l)/2. 

Notice that a and r commute, whereas aC, = (t and rf = £ a. Let (i. j) £ 
{0,1,..., n — l} 2 — {(0, 0)} be arbitrary. Then 

(a'r'Cf = [a‘ r 1 a - 1 ) = a i+j T i+j , 

hence 

(oV C) 2d = a {i+j)d T {i+j)d = (( ar ) i+i ) d , 

(<J*T J '0 2<i+ 1 = 0 -( i +J') d + i r ( i +J') d +J 

Clearly, (crV' 7 /) 2d+1 cannot be the identity permutation, so crV- 7 /' is of even 
order. Using the fact that, in the cyclic group {a, a 2 ,..., a™ -1 , a n = e} of order 
n, a k is of order n/gcd(fc,n), we find that the permutation i tV- 7 /) is of order 
2d, where d := n/ gcd(* + j. n). Therefore, every cycle of this permutation must 
have length that divides 2d. 

We claim that all cycles have length d or 2d. Accepting that for now, let 
us determine how many cycles have length d. A cycle that includes entry ( k,l ) 
has length d if (k,l) is a fixed point of {a l T^Q d . For this to hold we must have 
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d odd (otherwise there would be no fixed points because we have excluded the 
case i = j = 0 and (* + j)d /2 = lcm(* +j,n )/2 is not a multiple of n). Since 

{a i T i C) d = < j(’ + ^)( d_1 )/ 2+ *r^ +J '^ d_1 ^ 2+ - 7 'C, 

we must also have 


(&;, l) — ([Z + (i + j){d — l)/2 + j], [k + (i + j)(d — l)/2 + *]), (7) 

where d := n/ gcd(i + j, n) and, for simplicity, [r] := (r mod n ) £ {0,1,..., n — 
1}. For each k £ {0,1,..., n — 1}, there is a unique Z (namely, l := [fe + (i + 
j)(d — l)/2 + «]) such that (7) holds; indeed, 

[l + (i + j)(d — l)/2 + j] 

= [[& + (* + j){d — l)/2 + *] + (* + j)(d — l)/2 + j] 

= [k + (i + j)(d — l)/2 + i + (i + j){d — l)/2 + j] 

= [k + (i + j)d\ 

= [k + (i+j)(n/gcd{i+j,n))] 

= [k + lcm(* + j, n)] 

= k. 


This shows that there are n fixed points of {cr l T^Q d . Each cycle of length d of 
ij't^Q will account for d such fixed points, hence there are n/d such cycles. All 
remaining cycles will have length 2d, and so there are n(n — l)/(2 d) of these. 
The total number of cycles is therefore n{n + l)/(2 d). 

The other possibility is that d is even and all cycles have the same length, 2d, 
so there are n 2 /(2d) of them. Notice that d is a divisor of n, so the contribution 
to 

n— 1 n—1 

EE 2 "" 

i=0 j=0 


from odd d is 

rup{d) 2"("+ 1 )/( 2d ) 

d | n: d is odd 

and from even d is 

Y n<p(d) 2" 2 /( 2d ). 

d | n: d is even 


( 8 ) 

(9) 


The reason for the coefficient n<p(d) is that, if d | n, then the number of elements 
of the cyclic group {e, err, (err) 2 ,..., (err)” -1 } that are of order d is tp(d). And 
for a given ( i,j ) £ {0,1,..., n — l} 2 , there are n pairs (fc, l) £ {0,1,... ,n — l} 2 
such that [k + 1] = [i + j]. Putting (8) and (9) together, we obtain 


EE 2 

2=0 j= o 


Eij 


Y nv3(d)2 n( " +d_2Ld/2J)/(2d) , 

d | n 


( 10 ) 


which, together with (5) and (6), yields (3). 


□ 


6 


It remains to prove our claim that, for (*, j) S {0,1,..., n — l} 2 — {(0, 0)}, 
the permutation a l T J f cannot have any cycles whose length is a proper divisor 
of d := n/ gcd(i + j, n). Let c | d with 1 < c < d. We must show that (er 1 t^C,) c 
has no fixed points. We can argue as above with c in place of d. For (k,l) 
to be a fixed point of (a z T^f) c we must have ( i + j)c a multiple of n. But 
d := n/ gcd(i + j, n) is the smallest integer c such that (i + j)c is a multiple of 
n because ( i + j)n/ gcd(i + j , n) = lcm(i + j, n). 

Finally, we excluded the case n = 1 at the beginning of the proof, but we 
notice that the formula (3) gives a(l) = 2, which is correct. 


3 Rotation and reflection of rows and columns, 
and matrix transposition 


Let X n := {0, lfl 0,1 ’"'’" ^ be the set of n x n matrices of Os and Is, which 
has 2™ elements. Let /3(n) denote the number of orbits of the group action 
on X n by the group of order 8 n 2 generated by er (row rotation), r (column 
rotation), p (row reflection), 6 (column reflection), and f (matrix transposition). 
(Exceptions: If n = 2, the group is of order 8; if n = 1, the group is of order 1.) 

Informally, /3(n) is the number of (distinct) toroidal n x n binary arrays, 
allowing rotation and/or reflection of rows and/or columns as well as matrix 
transposition. 

Theorem 2. With b(n,n ) defined using (2), 0(n) is given by (4). 

Proof. Let us assume that n > 3. (We will treat the cases n = 1 and n = 2 
later.) By Polya’s enumeration theorem, 

^ n— 1 n— 1 

j3On) = —~ ^ y ( 2 Aij + 2 Bij + 2 Cij + 2 Dij + 2 Eij + 2 Fij + 2 Gij + 2 Hij ), 

8 n 2 ' 

2 — 0 J- =0 

where Aij (resp., Bij , C\j , D,j, Eij, Fij, Gij , Hjj ) is the number of cycles in the 
permutation ctV- 5 (resp., crV^p, crV J 0, a l T^pd, a a’d a pdf)-, 

here a rotates the rows (row 0 becomes row 1, row 1 becomes row 2, ..., row 
n — 1 becomes row 0), r rotates the columns, p reflects the rows (rows 0 and 
n—1 are interchanged, rows 1 and n — 2 are interchanged, ..., rows [_7T./2J — 1 
and n — |_n/2j are interchanged), 9 reflects the columns, and f transposes the 
matrix. The order of the group generated by a, r, p, 9, and f is 8n 2 , using the 
assumption that n > 3. 

We have already evaluated 

1 n—1n—1 

fl ( n ' n ) = ^EE 2 A «, 

i=0 j =0 


a(n) 


1 

27 ? 


y ^2(2 Aij + 2 Eij ), 

2=0 j —0 
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and 


b(n, n) = ^-2 ^2 E^' + + 2 Cij + 2° ij ), 


2—0 j — 0 


SO 


1 1 / 1 \ 1 n-ln-1 

P(n) = ^Hn,n)+-(a(n)--a(n,n)j+— I ^ ^(2 Fi *+2 G «+2 J?y ). (11) 

' ' 2—0 j— 0 


Let us begin with 


n—1n—1 

E E 2Hk 

i=0 j—0 


Here we are concerned with the permutations a l r°pOQ for (i,j) G {0,1,..., n — 
l} 2 . We will need some multiplication rules for the permutations er, t, p, 9 , and 
C, specifically 

err = rcr, (70 = 9(7, rp = pr , p9 = 9p, ap = pcr~ 1 , t 9 = dr^ 1 , 

and 

crC = C t j < = O, pC = C0, 0 C = Cp- 

It follows that (with t~ 1 := (r -1 ) 2 ) 

a l r^ pOC, = a l r^C,6p = (T^a-'Op = (9r~ l a-’ p = C,8pT~ l a~\ 


and hence 

(aV-VC) 2 = (c7Vp6>C)(C0pT _ V-O = a l ~ 3 T~ l+j = (ot-- 1 ) 4-5 ' = (o-^tP+L 

(12) 

In particular, if i 6 (0,1,..., n— 1}, then the permutation a l T l p9( is of order 
2. Furthermore, under this permutation, the entry in position (k,l) moves to 
position (n — 1 — [Z + i], n — 1 — [fc + «]), where, as before, [r] := (r mod n) G 
(0,1,..., n — 1}. Thus, (k, Z) is a fixed point if and only if 

(k, l) = (n — 1 — [l + *], n — 1 — [k + i]). (13) 

For each k G {0, 1, ..., n — 1} there is a unique l G {0,1,..., n — 1} (namely 
l := n — 1 — [k + *]) such that (13) holds; indeed, 

n — 1 — [l +i\ = n — 1 — [n — 1 — [k + i]+i]=n — 1 — [n — 1 — (k + i) + i] 

= n — 1 — [n — 1 — k\ = n — 1 — (ro — 1 — k) = k. 

Thus, a l T l p9( with i G {0,1,..., n — 1} is of order 2 and has exactly n fixed 
points, hence (’ 2 l ) transpositions. This implies that Ha = n(n + l)/2 for such i. 

Now we let (i,j) G {0,1,... ,n — l} 2 be arbitrary but with i ^ j. Let us 
generalize (12) to 

(o-Vp0C) 2d = a (l - j)d T ( - l+j)d = ((aT~ 1 ) i - j ) d = ((a~ 1 T)- i+j ) d , 
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(oVp6»C) 2d+1 = a^- j)d+i T ( - l+j)d+j p6(. 

The proof proceeds much like the proof of Theorem 1. Specifically, a l r 3 pOC, is 
of order 2d, where d := n/ gcd(|* — j\,n). All cycles have length d or 2d. In 
fact, if d is odd, there are n/d cycles of length d and n{n — 1)/(2c?) cycles of 
length 2d. If d is even, there are n 2 /(2d) cycles, all of length 2d. And for a 
given [i, j) £ {0,1,..., n — l} 2 , there are n pairs ( k , l ) £ {0,1,..., n — l} 2 such 
that [k — l] = [| i — j |]. We arrive at the conclusion that 

n— 1 7i—l n—1 71—1 

J2J2 2Hij = J2J2 2Eij - ( l4 ) 

i=0 j=0 i=0 j —0 

Next we evaluate 

71—1 71—1 71—1 71—1 

EE 2Fy = EE 2Gy ’ ( 15 ) 

i=0 j =0 i=0 j —0 

where the equality holds by symmetry. We consider the permutations a l r J pC, 
for (i,j) £ {0,1 l} 2 . From the multiplication rules, it follows that 

crV- 7 pQ = 


and hence 


(crV- 7 p() 2 = (crV- 7 p(){(9t V j ) 
which implies 


(jV J p0T = a 1 ^t 1+ ^ pO = 9pr 1 a I+J , 

(16) 


(o-VVC) 4 = (<r i- V +J ' p0) {Opr-*--i(j- i+j ) = e. 


So the permutation <j 1 t 3 pC, is of order 4. The entry in position (k,l) moves 
to position ([l + j],n — 1 — [k + i\) under this permutation. Thus, (k,l) £ 
(0,1,..., n — l} 2 is a fixed point of <j 1 t 3 pC, if and only if 


(M) = {[l+j],n- 1 — [k + i]). 

There is a solution (fc, l) if and only if there exists l £ {0,1,..., n— 1} such that, 
with k := [/ + j], we have n — 1 — [k + i] = l or, equivalently, 

[l + * + j) = n — 1 — l. (17) 


When i + j < n — 1, (17) is equivalent to 

l + i + j = n— 1 — l or l + i + j — n = n — 1 — l 


or to 

l = (n — 1 — i — j)/2 or l = (2n — 1 — i — j)/2. 


If n is odd and i+j is odd, then there is one fixed point, (fc, l) = ([(2 n —1 —i + 
j)/ 2], [(2n — 1 — i — j)/2]). If n is odd and i + j is even, then there is one fixed 
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point, (k, l) = ([(n — 1 — i + j)/ 2], [(n — 1 — i — j)/ 2]). If n is even and i + j is 
odd, then there are two fixed points, namely 

(M) = ([(«- 1 ~i+j)/ 2], [(n - 1 - i-j)/ 2 ]), 

(M) = ([(2n — 1 — * +i)/2], [(2n — 1 - i - j)/2\). 

Finally, if ro is even and i+j is even, then there is no fixed point. 

When i + j > n, (17) is equivalent to 

l + i + j ~ n = n — 1 — l or l + i + j — 2n = n — 1 — l 

or to 

l = (2n — 1 — i — j)/2 or l = (3n — 1 — i — j)/2. 

If n is odd and i+j is odd, then there is one fixed point, (k, l) = ([(2 n — 1 — i + 
j)/ 2], [(2n—1— i— j)/2]). If n is odd and i+j is even, then there is one fixed point, 
(fc, l) = ([(n-l-i+j)/2], [(3n—1—i—j)/2]) = ([(n-l-i+j)/2], [{n-l-i-j)/2}). 
If n is even and i + j is odd, then there are two fixed points, namely 

(M) = ([(2 n - 1 -i + j)/ 2], [(2n — 1 - i - j)/2]), 

(*, 0 = ([(^ — 1 — * + j)/2], [(^ — 1 — * — j')/2]). 

Finally, if n is even and i+j is even, then there is no fixed point. Notice that 
the results are the same for i + j > n as for i + j < n — 1. 

Using (16), under the permutation (crV-lpC) 2 , the entry in position (k,l) 
moves to position (n — 1 — [k + i — j], n — 1 — [l + i + j]). Thus, (k,l) £ 
{0,1,..., n — l} 2 is a fixed point of (cr'r^pC) 2 if an( i on ly if 

(k,l) = (n - 1 - [k + i - j],n - 1 - [l + i+j]). 

A necessary and sufficient condition on (k,l) is (17) together with [k + i — j] = 
n — 1 — k. Solutions have l as before. On the other hand, k must satisfy 

k + i — j — n = n— 1 — k, k + i — j = n— 1 — k, or k + i — j + n = n— 1 — k, 

or equivalently, 

k = [(n — 1 — i + j)/2] or k = [(2n — 1 — i + j)/2\. 

If n is odd, the only fixed points of pC,) 2 are those already shown to be 
fixed points of cr l T J pC,. If n is even and i + j is odd, there are two fixed points 
of (aV-ipC) 2 that are not fixed points of ++p(j. namely 

(M) = ([(n - 1 - i + j)/2],[(2n - 1 - i - j)/2]), 

(■M) = ([(2 n - 1 -i + j)/ 2], [(n- 1 -i -j)/2]). 

Finally, there are no fixed points when n is even and i+j is even. 
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Consequently, if n is odd, then the permutation <j 1 t 3 p£, which is of order 4, 
has only one fixed point. Therefore, it has one cycle of length 1 and ( n 2 — l)/4 
cycles of length 4. Thus, 


EE 2Fii 

2=0 j = o 


= n 2 2(” 2 + 3 )/ 4 . 


For even n, if i+j is odd, then the permutation a 1 t° pC, has two cycles of length 
1 and one cycle of length 2, and the remaining cycles are of length 4. If i + j is 
even, then all cycles of the permutation a l T 3 pC, are of length 4, hence there are 
tt 2 /4 of them. Thus, 


E E 2Fij 

2=0 j —o 


I n 2 2 (n 2 - 4 )/ 4 + 3 I n 2 2 n 2 /4 = 5n 2 2 ™ 2 /4-l_ 

2 2 


These results, together with (3), (10), (11), (14), and (15), yield (4). 

Finally, recall that we have assumed that n > 3. We notice that the formula 
(4) gives /3(1) = 2 and /?(2) = 6, which are correct, as we can see by direct 
enumeration. □ 


In the derivation of (2) in [1], the proof requires to, n > 3 because the group 
D m x D n used in the application of Polya’s enumeration theorem ( D m being 
the dihedral group of order 2 m), is incorrect if to or n is 1 or 2. If m = 2, row 
rotation and row reflection are the same, so the latter is redundant. Thus, Di 
should be replaced by C2, the cyclic group of order 2. The reason (2) is still 
valid is that b 1 (2,n) = 62(2,71) and 63(2,71) = 64(2,71), as is easily verified. If 
m = 1, again row reflection is redundant, so D\ should be replaced by C\. Here 
(2) remains valid because 61 (1,71) = 62(1,71) and 63(1,71) = 64(1,71). A similar 
remark applies to n = 2 and n = 1, except that here 61 (to, 2) = 63(771,2), 
62(771, 2) = 64(771, 2), 61 (m, 1) = 63(771,1), and 62(771,1) = 64(777,1). 
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